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Abstract 



Recent theoretical advances in the study of heavy ion fusion reactions below 
the Coulomb barrier are reviewed. Particular emphasis is given to new ways 
of analyzing data, such as studying barrier distributions; new approaches to 
channel coupling, such as the path integral and Green function formalisms; 
and alternative methods to describe nuclear structure effects, such as those 
using the Interacting Boson Model. The roles of nucleon transfer, asymmetry 
effects, higher-order couplings, and shape-phase transitions are elucidated. 
The current status of the fusion of unstable nuclei and very massive systems 
are briefly discussed. 
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I. INTRODUCTION 



Quantum tunneling in systems with many degrees of freedom is one of the fundamental 
problems in physics and chemistry (Hanggi et ai, 1990; Tsukada et ai, 1993). One example 
of tunneling phenomena in nuclear physics is the fusion of two nuclei at very low energies. 
These reactions are not only of central importance for stellar energy production and nucle- 
osynthesis, but they also provide new insights into reaction dynamics and nuclear structure. 
Until about fifteen years ago, low energy fusion reactions were analyzed in terms of a simple 
model, where one starts with a local, one-dimensional real potential barrier formed by the 
attractive nuclear and repulsive Coulomb interactions and assumes that the absorption into 
the fusion channel takes place at the region inside the barrier after the quantum tunnel- 
ing. The shape, location, and the height of this potential were described in terms of few 
parameters which were varied to fit the measured cross sections. The systematics of poten- 
tials obtained in this way were discussed by Vaz et ai (1981). A number of experiments 
performed in the early 80's showed that the subbarrier fusion cross sections for intermediate- 
mass systems are much larger than those expected from such a simple picture (Beckerman, 
1988; Vandenbosch, 1992). The inadequacy of the one-dimensional model for subbarrier 
fusion was explicitly demonstrated by inverting the experimental data to directly obtain the 
effective one-dimensional fusion barrier under the constraint that it is energy independent 
(Balantekin et ai, 1983). 

In recent years there has been a lot of experimental effort in measuring fusion cross 
sections and moments of compound nucleus angular momentum distributions. A complete 
compilation of the recent data is beyond the scope of this review. Several recent reviews 
(Stefanini et ai, 1994; Reisdorf, 1994) present an excellent overview of the current exper- 
imental situation. This article reviews theoretical developments of the last decade in our 
understanding the multidimensional quantum tunneling nature of subbarrier fusion. We 
should emphasize that the large volume of new data and extensive theoretical work does 
not allow us to provide an exhaustive set of references. The selection we made does not 
imply that omitted references are any less important than the ones we chose to highlight in 
discussing different aspects of fusion phenomena. 

The natural language to study fusion reactions below the Coulomb barrier is the coupled- 
channels formalism. In the last decade coupled channels analysis of the data and the re- 
alization of the connection between energy derivatives of the cross section and the barrier 
distributions (Rowley et ai, 1991) motivated accumulation of very high precision data. 

When the enhancement of the cross section below the barrier was first observed, many 
authors pointed out that it is not easy to identify the underlying physical mechanism (Brink 
et ai, 1983b; Krappe et ai, 1983). Any coupling introduced between translational motion 
and internal degrees of freedom enhances the cross section. The recent high-precision data 
helped resolve some of these ambiguities by studying barrier distributions and mean angular 
momenta as well. The quality of the existing data now makes it possible to quantitatively 
explore a number of theoretical issues in the quantum tunneling aspects of subbarrier fusion, 
such as effects of anharmonic and non-linear interaction terms. 

In the next section, after briefly reviewing observables accessible in heavy-ion fusion 
reactions, and the reasons why a one-dimensional description fails, we present the motivation 
for studying barrier distributions. In Section |T| we discuss the standard coupled-channels 
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formalism, and alternative approaches such as the path integral formalism and Green's 
function approaches along with their various limiting cases. Section [V has a dual purpose: 
it covers recent attempts to describe nuclear structure effects using the Interacting Boson 
Model while illustrating the technical details of these alternative approaches discussed in 
Section [ITj|. In Section [V], we survey a representative sample of recent high-quality data 
with a focus on new physics insights. We conclude in Section |Vj with a brief discussion of 
the open problems and outlook for the near future. 



II. ONE-DIMENSIONAL MODEL FOR FUSION 

A. Experimental Observables 

In the study of fusion reactions below the Coulomb barrier the experimental observables 
are the cross section 

oo 

a(£)=J>(£), (1) 



and the average angular momenta 



The partial-wave cross sections are given by 

a e (E) = ^(2£+l)T e (E), (3) 

where T^{E) is the quantum- mechanical transmission probability through the potential bar- 
rier and \i is the reduced mass of the projectile and target system. 

Fusion cross sections at low energies are measured by detecting evaporation residues or 
fission products from compound nucleus formation. Evaporation residues can be detected 
directly by measuring the difference in velocities between them and the beam-like ions. 
Velocity selection can be achieved either by electrostatic deflectors or by velocity filters. 
Alternative techniques include detecting either direct or delayed X-rays and gamma rays. A 
review of different experimental techniques to measure the fusion cross sections is given by 
Beckerman (1988). 

Several techniques were developed for measuring moments of the angular momentum 
distributions. The advent of detector arrays enabled measurement of full gamma-ray mul- 
tiplicities (Fischer et ai, 1986; Halbert et ai, 1989), which, under some assumptions, can 
be converted to oy distributions. It is also possible to measure relative populations of the 
ground and isomeric states in the evaporation residues to deduce the spin distribution in the 
compound nucleus (Stokstad et ai, 1989; DiGregorio et ai, 1990). Finally, the anisotropy 
of the fission fragment angular distribution can be related to the second moment of the spin 
distribution (Back et ai, 1985; Vandenbosch et ai, 1986a; Vandenbosch et ai, 1986b). 

It is worthwhile to emphasize that moments of angular momenta, unlike the fusion cross 
section itself, are not directly measurable quantities. One needs to make a number of 
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assumptions to convert gamma-ray multiplicities or isomeric state populations to average 
angular momenta. An elaboration of these assumptions along with a thorough discussion of 
the experimental techniques is given by Vandenbosch (1992). 



B. One-Dimensional Barrier Penetration 



The total potential between the target and projectile nuclei for the £-th partial wave is 
given by 



VAr) = V N (r) + Voir) + 



h 2 i(£ + 1) 
2/ir 2 



Voir) + 



2/ir s 



(4) 



where Vjv and Vc are the nuclear and Coulomb potentials, respectively. The i = barrier is 
referred to as the Coulomb barrier. The barriers obtained in Eq. (|j) is illustrated in Figure 
[I] for several i- values. 

For a one-dimensional barrier transmission probabilities can be evaluated numerically - 
either exactly or using a uniform WKB approximation, valid for energies both above and 
below the barrier (Brink and Smilansky, 1983a; Brink, 1985a) : 



T e (E) = [l + exp(2S e (E))}- 1 
where the WKB penetration integral is 



Si(E) 




2/X f r M 



dr 



'it 



V (r) 



h 2 



2/ir 2 



— E 



1/2 



(5) 



(6) 



In this equation and r 2 i are the classical turning points for the £-th partial wave potential 
barrier. 

If we assume that the potential barrier can be replaced by a parabola 



V ir) = V B o-y 2 n 2 ir-r f 



(7) 



where Vbo is the height and Q is a measure of the curvature of the s-wave potential barrier, 
the transmission probability can be calculated to be (Hill and Wheeler, 1953) 



To(E) 



1 + exp 



2tt , 



(8) 



In the nuclear physics literature Eq. (|j) is known as the Hill- Wheeler formula. Especially 
at energies well below the barrier there are significant deviations from this formula as the 
parabolic approximation no longer holds. 



C. Barrier Distributions 

An alternative way of plotting the total cross section data is to look at the second 
energy derivative of the quantity Ecr, sometimes called the distribution of the barriers. To 
elaborate on the physical significance of this quantity let us consider penetration probabilities 
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for different partial waves in the case of a one- dimensional system (coupling to an internal 
system is neglected), Eq. (g). Under certain conditions, to be elaborated in the next 
subsection, we can approximate the I dependence of the transmission probability at a given 
energy by simply shifting the energy (Balantekin et al, 1983; Balantekin and Reimer, 1986): 



E- 



£(£ + l)h 2 
2/jR 2 (E) 



(9) 



where fiR 2 (E) characterizes an effective moment of inertia. R(E) was found to be a slowly 
varying function of energy as depicted in Figure |!|. Consequently, in many applications 
R(E) is replaced by tq, the position of the s-wave barrier, in Eq. (|9|). If many values of 
I are important in the sum over partial-wave transmission probabilities in Eq. (1), we can 
approximate that sum with an integral over £, and, using Eq. (|S|) obtain (Balantekin et al, 
1983; Balantekin and Reimer, 1986) 

Ea(E) = nR 2 (E) [ E dE'T (E'). (10) 



It was found that Eq. ( |10|) represents the experimental data for the total fusion cross 
section rather well (Balantekin and Reimer, 1986; Dasgupta et al, 1991; Balantekin et 
al, 1996). Differentiating Eq. ( |T0D twice (Rowley et al, 1991) one finds that the energy 
derivative of the s-wave transmission probability is approximately proportional to the second 
energy derivative of the quantity Ea up to corrections coming from the energy dependence 
of R(E): 

dTn(E) 1 cP .„ ,„,dR, 

Since R(E) is a slowly varying function of energy, the first term in Eq. ( |TTD can be used 
to approximate the first derivative of the s-wave tunneling probability. For a completely 
classical system, T is unity above the barrier and zero below; hence the quantity dT (E) / dE 
will be a delta function peaked when E is equal to the barrier height, as shown in Figure 
[3]. Quantum mechanically this sharp peak is broadened as the transmission probability 
smoothly changes from zero at energies far below the barrier to unity at energies far above 
the barrier (see Figure ||). Rowley et al. (1991) suggested that if many channels are coupled 
to the translational motion, the quantity dT Q (E) jdE is further broadened and can be taken 
to represent the "distributions of the barriers" due to the coupling to the extra degrees of 
freedom as depicted in Figure |] for a two-channel case. 

Obviously one needs rather high precision data to be able to numerically calculate the 
second derivative of the excitation function. Such high-quality data recently became avail- 
able for a number of systems. As a specific example of the quality of recent data, Figure |5| 
shows the measured cross section and the extracted associated barrier distribution for the 
16 + 154 Sm system (Wei et al, 1991; Leigh et al, 1995). 

Eqs. (H) and (|T0| ) (with the substitution R(E) — > r ) can also be used to obtain a direct 
connection between the fusion cross section and the angular momentum distribution (Sahm 
et al, 1985; Ackermann, 1995) : 
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1 (d(E'a(E')Y 



dE' 



(12) 



where 



E' = E- (h 2 /2/ir 2 )£(£ + l). 



(13) 



For energies well above the barrier, one can use the parabolic approximation of Eq. (^) 
for the potential barrier. Further approximating R(E) by r and inserting the penetration 
probability for the parabolic barrier, Eq. (^j), in Eq. fllOD one obtains an approximate 
expression for the cross section (Wong, 1973) 



a{E) 



Tittr 2 
2E 



log 



1 + exp 



2tt 

m 



(E — Vi 



BO) 



(14) 



In the classical limit, where Q — > or E ^> Vbo, Eq. ( |T4|) reduces to the standard geometrical 
result 



(15) 



D. Energy Dependence of the Effective Radius 

If one sets R(E) = in Eq. (H) for approximating the £-wave penetrability by the 
s-wave penetrability at a shifted energy, one gets only the leading term in what is actually 
an infinite series expansion in A = £[£ + 1). The next term in this expansion can easily be 
calculated. Let r e denote the position of the peak of the £-wave barrier which satisfies 



dV t (r) 



and 



Or 



d 2 VAr) 







(16) 



dr 2 



< 



(17) 



r=r e 



then the height of the barrier is given by Vbi = VJ(r^). We make the ansatz that the barrier 
position can be written as an infinite series, 



re = r + ciA + c 2 A 2 + 



(18) 



where the c« are constants. Expanding all functions in Eq. ( |i~6"D consistently in powers of A, 
we find that the first coefficient is 



h 2 



ci 



(19) 



where a is the curvature of the s-wave barrier 
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a = — 



d 2 M 
Or 2 



Substituting the leading order correction in the barrier position re into Eq. 
to second order in A the £-wave barrier height is given by 

v =v | h " A | n ' A2 

B B0 2jj,rQ 2^i 2 ar% ' 



(20) 



), we find that 



(21) 



Therefore, an improved approximation for the ^-dependence in the penetrability is given by 

i h 2 A h A A 2 \ 



T e (E) « T [E 



2fir 2 



2 / u 2 aro / 



(22) 



Comparing Eq. (21) with Eq. (6), we find that the energy-dependent effective radius can be 
expressed as (Balantekin et al, 1996) 



R 2 {E) 



4 J E dE'T (E')(E-E') 



tfdE>T (E>) 



(23) 



This expression is useful in assessing the utility of the second derivative of the quantity Eo 
as the distribution of barriers as it gives an estimate of the terms neglected in Eq. (|TT|) . If 
we rewrite Eq. ( |TTD including previously neglected terms 



d 2 



dT (E) _ 1 

dE 7rR 2 (E)dE 2 



(ME)) 



Ea(E) d 2 
ttR^(E) dE 2 



(R 2 (E) 



2T (E) d 
R 2 (E) dE 



(R 2 (E)) , (24) 



we see that a strong energy dependence of R(E) would not only provide an overall mul- 
tiplicative factor between the experimental observable d 2 (Ea)/dE 2 and the true barrier 
distribution (i.e., dT /dE), but may also shift the position of various peaks in it and change 
the weights of these peaks. Eq. ( p3|) can be used to illustrate that such corrections are 
indeed small (Balantekin et al, 1996). This is a useful consistency check even when channel 
coupling effects yield a number of eigen-barriers (cf. Section |llllj| ) as Eq. QTTj) still needs 
to be satisfied for each one-dimensional eigen-barrier to be able to interpret experimentally 
determined d 2 (Ea)/dE 2 as the distribution of barriers. 



E. Inversion of the Data 

Using Abelian integrals (Cole and Good, 1978), one can show that for energies below 
the barrier 



E 



y/E'-E 2\hf 



V (r) + 



h 2 £(£ + 1) 
2{ir 2 



E 



(25) 



where Vbi is the height of the £-wave potential and Vq{t) is the s-wave barrier. The energy 
derivative of the left-hand side of Eq. (53) can be integrated by parts to yield, 
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y B t dE , dS e (E')/dE' 

E y/E'-E 




(26) 



which is used to find the barrier thickness (Balantekin et al, 1983). Using Eqs. (BT) and (^j), 
one can relate the WKB penetration integral to the experimentally measured cross section 



as 



£„(£) = -log 



_d_ 

dE 



EajE) 
txR 2 (E) 



(27) 



Thus, if R(E) is specified, the thickness of the barrier at a given energy is completely 
determined from the experimental data using Eqs. ( [ZED an d (p7p. 

The potentials resulting from the analysis of Balantekin et al. (1983) for six systems are 
shown in Figure || For comparison, point Coulomb and the phenomenologically determined 
potential of Krappe et al. (1979) are also shown. It should be emphasized that this inversion 
method assumes the existence of a single potential barrier. The thickness functions, t = 
^20 — r io ; m Figure ^| especially for the heavier systems are inconsistent with the assumption 
of a single-valued one-dimensional local potential, clearly indicating the need for coupling 
to other degrees of freedom. This result was confirmed by the systematic study of Inui and 
Koonin (1984). 



III. MULTIDIMENSIONAL QUANTUM TUNNELING IN NUCLEAR PHYSICS 

A. Coupled-Channels Formalism 



A standard theoretical approach to study the effect of nuclear intrinsic degrees of freedom 
on the fusion cross section is to numerically solve coupled-channels equations that determine 
the wave functions of the relative motion. For example, if one is interested in the effect of 
the excitation of the ground state K = + rotational band of the target nucleus, then 
each channel can be denoted by a set of indices (I,£), where I and I are the angular 
momentum of the rotational excited states of the target nucleus and that of the relative 
motion, respectively. The coupled-channels equations then read 

h 2 d 2 h 2 



2\i dr 2 2/ir 2 



(2£ 1 + l)(2/ 1 + l)(2/ 2 + l) rl 1 



2A + 1 



9 ihJ] 



x < A/200IA0 >< £iA00|40 > { jj >fh,i 2 (r) = 0. 



(28) 



When the coupled-channels formalism is used to study direct reactions the optical potential 
V(r) contains an imaginary part in order to take into account the effect of intrinsic degrees 
of freedom other than the rotational excitation on the scattering. The multi-polarity of the 
intrinsic excitation is A = 2 if we restrict to quadrupole deformations. The coupling form 
factor F\(r) consists of Coulomb and nuclear parts, 



9 



F x (r) = F c (r)+F N (r). 



(29) 



If the coupling is restricted to the quadrupole deformations, the Coulomb part is 
F c (r) = -L=/3 Zl Z 2 e^ (r > R c ) 



3 --PZiZ 2 e 2 ^ (r<R c ) (30) 



V20^ R 3 C 



and the nuclear part is 



F N {r) = -\l^mv^ ■ (31) 



In most calculations the scale parameters R c and Ry are taken to be 1.2 A 1 ^. In the standard 
coupled-channels calculations, one solves Eq. ( f2~8~D by imposing regular boundary condition 
at the origin. In contrast, in the study of heavy- ion fusion reactions, one takes the potential 



to be real and often solves Eq. (28) by imposing the incoming wave boundary condition at 
some point inside the potential barrier to obtain S-matrices. The fusion cross section is then 
obtained based on the unitarity relation as 



7T 

F 2 



a = ^£(2£ + l){l-£|S a (£)| 2 }, (32) 



where the index a designates different scattering channels, which have been explicitly dealt 
with in the coupled-channels calculation. An advantage of the coupled-channels method is 
that one can try to consistently analyze heavy-ion fusion reactions, and other scattering pro- 
cesses such as elastic scattering. Coupled-channels calculations for a number of systems exist 
in the literature (Dasso et al, 1983; Lindsay and Rowley, 1984; Esbensen and Landowne, 
1987; Thompson et al, 1985; Stefanini et al, 1990). 

One can significantly reduce the dimension of the coupled-channels calculations by ig- 
noring the change of the centrifugal potential barrier due to the finite multi-polarity of 
the nuclear intrinsic excitation (Takigawa and Ikeda, 1986; Esbensen et al, 1987). This is 
called the no-Coriolis approximation, rotating frame approximation, or iso-centrifugal ap- 
proximation (Gomez-Camacho and Johnson, 1988). A path integral approach to no-Coriolis 
approximation was given by Hagino et al. (1995a). 

If one further ignores the finite excitation energy of nuclear intrinsic motion, then one 
can completely decouple the coupled-channels equations into a set of single eigen-channel 
problems. These two approximations significantly simplify the numerical calculations and 
also give a clear physical understanding of the effect of channel-coupling in terms of the 
distribution of potential barriers. 

B. Simplified Coupled-Channels Models 

Under certain assumptions it is possible to significantly simplify the coupled-channels 
equations described in the previous section. We take the Hamiltonian to be 
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H = H k + V (r) + H (0 + H int (r, f) 



(33) 



with the term H k representing the kinetic energy 

h 2 , 

H k = V 2 

2/i 



(34) 



where r is the relative coordinate of the colliding nuclei and £ represents any internal degrees 
of freedom of the target or the projectile. In this equation Vo(r) is the bare potential and 
the term H (£) represents the internal structure of the target or the projectile nucleus. 
Introducing the eigenstates of H (£) 



H \n) = e n \n), 
and expanding the radial wave function as 



(35) 



(36) 



the time-independent Schroedinger equation is reduced to a set of coupled equations for the 
relative motion wave functions Xn, 



+ v ' ir) 



E 



Xn(r) = - YJtnSnm + (n \ H lnt (r, £) | TTl) ] Xm (0 
m 

These equations are solved under the incoming wave boundary conditions: 

5 n0 exp(-ik n r) + x PfR n exp(ik n r), r -> +00 



(37) 



Xn(r) 



jrT n exp(-ik n r), 



(38) 



where h 2 k 2 l /2fi = E — e n . The internal degrees of freedom are taken to be initially in their 
ground state labeled by n = and the associated ground state energy is set to be zero, 
eo = 0. In Eq. (|38D the reflection and transmission coefficients in each channel are denoted 
by R n and T n , respectively. 

Several groups (Dasso et al, 1983; Broglia et al, 1983a; Jacobs and Smilansky, 1983) 
studied various simplifying limits to emphasize salient physical features. Here we summarize 
the approach of Dasso et al. (1983). They assumed that the coupling interaction factors into 
a relative part, F(r), and an intrinsic part, 



M nm = e n 5 nm + (n\H int (r, £)|ra) = e n 5 nm + F(r)(n\G(g)\m), 



(39) 



and that the form factor F(r) is a constant (taken by Dasso et al. (1983) to be the value 
of F(r) at the barrier position). Under these approximations the coupled-channel equations 
decouple to give 



2/i dr 2 



+ VAr) + \ n -E 



Y. U nmXm(r) = 0, 



(40) 
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where U nm is the unitary matrix which diagonalizes the coupling matrix M nm to give a set 
of eigenvalues A n . Eq. (|40|) indicates that the effect of the coupling is to replace the original 
barrier by a set of eigen-barriers Vf(r) + A n . The transmission probability calculated in the 
incoming-wave boundary conditions is given by 

T e (E)=J2\U n0 \ 2 Te(E-\ n ), (41) 

n 

where Te(E — A n ) is the transmission probabilities calculated at shifted energies E — A n . 

Even though the const ant- coupling approach would overpredict the transmission prob- 
ability, it can nevertheless be used to get a qualitative understanding of the dependence of 
the fusion cross section on various physical quantities. For example, one can study coupling 
to an harmonic mode with a finite Q-value using the model Hamiltonian 

M = tt 2 /2D + l -Ce + F £, (42) 

yielding 



M mn = -nQS mn + F(y/nS n>m+1 + y/n + l£ n>m _i), (43) 



where — Q = hyC/D is the excitation energy and F = F ^\Q\/2C measures the total 
strength of the coupling and the eigenvalues are 

X n = n\Q\-F 2 /\Q\. (44) 

The total transmission probability can easily be calculated to be 

oo 

T CC (E) = Y,(F 2n /Q 2n n\)exp(-F 2 /Q 2 )T(E-n\Q\+F 2 /\Q\). (45) 

n=0 

In the special case of a two-channel problem the matrix 

F 

M=\ I (46) 



Q 



has eigenvalues 



A ± = -(-Q±VQ 2 + 4F2), (47) 



with the corresponding weight factors 

2F 2 

u2± = 4F* + Q*tQ^ + Q^ (48) 

Note that for F/\Q\ < 1 the lowest effective barrier carries the largest weight for negative 
Q, while the situation is reversed for positive Q (cf. Section |V U| ). 

Tanimura et al. (1985) pointed out that the constant coupling approximation can over- 
estimate low-energy fusion rates where the coupling effects are strong and the associated 
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form factors are rapidly varying. Dasso and Landowne (1987a) extended their model for 
strong coupling cases. They diagonalized the coupled intrinsic system in the barrier region 
to obtain the eigenstates | j(r)) and the eigenvalues A 7 (r) as functions of r to obtain the 
total transmission probability 

T SCC (E) = £ | ( 7 (J2) | 0}\%(E, V(r) + A 7 (r)), (49) 

7 

where T y (E, V{r) + A 7 (r)) is the penetration probability for the potential V(r) + A 7 (r). The 
weighing factors are fixed at a chosen radius R, which might be the position of the unper- 
turbed barrier or the average position of the eigen-barriers. Two coupled-channels codes, 
simplified in this manner, are available in the literature: CCFUS (Dasso and Landowne, 
1987b) and CCDEF (Fernandez-Niello et ai, 1989). Both of them treat the vibrational cou- 
pling in the constant coupling approximation. The latter takes into account projectile and 
target deformations within the sudden approximation and treats coupling to the transfer 
channels with a constant form factor. 

In CCFUS the basis states included are the ground state | 0), the quadrupole one-phonon 
state &r, | 0), the octupole one-phonon state b 3 | 0), and the product two-phonon state 6363 | 0). 
The resulting matrix to be diagonalized to yield the eigen-channels is 

M ',:";[ ',; " . (50) 



/ F 2 {r) F 3 (r) \ 

F 2 (r) e 2 F 3 (r) 

F 3 (r) e 3 F 2 (r) 

\ F 3 (r) F 2 (r) e 2 + e 3 J 



In CCFUS the double phonon states {b\f | 0), and {b\f | 0) are not included for mathematical 
simplification. The eigenvalues of the matrix in Eq. (|50"D can be written as the sums of the 
eigenvalues of the two 2x2 matrices which represent the coupling of the single one-phonon 
states, i.e., 



M2 =Um h )• (51) 

and 

If the matrix of Eq. (|50|) also included the double-phonon states (6 2 ) 2 |0), and (b\) 2 | 0) 
such a simple relationship between eigenvalues of the 4x4 and 2x2 matrices would not 
be possible. CCFUS provides two options; the matrices of Eq. (|5T|) and (|52|) are either 



diagonalized by replacing the form factors F 2 (r) and F 3 (r) with their values at the location 
of the bare potential barrier, or diagonalized for all values of r. In the latter case only the 
weight factors, but not the eigenvalues are evaluated at the position of the bare barrier. The 
transmission coefficients for each eigen-barrier are calculated using the Wong formula, Eq. 

(0)- 

The approach of CCFUS can be generalized to incorporate n different phonons by in- 
cluding the ground state, n one-phonon states, n(n — l)/2 two-phonon states (i.e., not those 
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states where the same phonon appears more than once), n(n — l)(n — 2)/3! three-phonon 
states, etc. to obtain a total number of 



1 + n + 



n(n — 1) n(n — l)(n — 2) 



+ 



n(n — l)(n — 2)(n — 3) 
4! 



+ 



2 



(53) 



2! 3! 



states. Since no phonon appears more than once, the eigenvalues of the resulting 2 n x 2 n 
matrix can be written as appropriate combinations of the eigenvalues of n 2 x 2 matrices. 
Even though the approach of CCFUS provides an elegant mathematical solution to the 
matrix diagonalization problem it ignores all the states where the same phonon appears 
more than once, e.g., the double-phonon states. Dasgupta et al. (1992) and Kruppa et al. 
(1993) pointed out that in some cases the coupling of a state like (fej) 2 | 0) to the ground 
state can be stronger than the coupling of a state like b\bj | 0} , i ^ j. Kruppa et al. (1993) 
considered quadrupole and octupole phonons, but included their double-phonon states as 
well diagonalizing the resulting 6x6 matrix. Dasgupta et al. (1992) excluded all multiple- 
phonon states, so for n different types of phonons they numerically diagonalized an (n + 
1) x (n + 1) matrix instead of analytically diagonalizing an 2 n x 2 n matrix. The resulting 
simplified coupled channels code is named CCMOD (Dasgupta et al, 1992). In CCMOD the 
weight factors are calculated at the position of the bare barrier, but the energy dependence 
of R(E) in Eq. (||) is taken into account using the prescription of Rowley et al. (1989). 
Eigen-channel cross sections are again calculated using the Wong's formula. 

Many experimentalists use these simplified coupled channel codes. When using them it 
is important to remember the approximations discussed in the previous paragraphs. Some 
of these approximations (constant coupling, Wong's formula) lead to an overestimate of the 
cross section. Some of the ignored couplings (e.g., the double phonon states) may be very 
important for the dynamics of the analyzed system. We therefore recommend using these 
codes only for a qualitative understanding of the data and strongly encourage authors to 
use full coupled channel codes for any quantitative description. 

Finally one should point out that in the limit in which the intrinsic energies are small 
compared to the coupling interaction one can approximate Eq. ( |39| ) 



In this case it is not necessary to require the coupling to be constant. The transformation 
amplitude between the ground state and the eigen-channels labeled by £ is the ground state 
wave function yielding the transmission probability 



where T[E, V(r)+F(r)G(^)} is the transmission probability for the potential V(r)+F(r)G(£) 
calculated at energy E. 



An alternative formulation of the multidimensional quantum tunneling is given by the 
path integral formalism (Balantekin and Takigawa, 1985). For the Hamiltonian given in Eq. 




(54) 




(55) 



C. Path Integral Approach 
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(|33|) the propagator to go from an initial state characterized by relative radial coordinate 
(the magnitude of r) and internal quantum numbers rij to a final state characterized by 
the radial position rf and the internal quantum numbers n/ may be written as 



K{r f ,n f ,T;r hni ,0) = / V [r(t)) e^ r ^W nfni (r(t), T) 



(56) 



where S(r,T) is the action for the translational motion and W nfni is the propagator for the 
internal system along a given path of the translational motion: 



W nfnt (r,T) = 

Ui n t satisfies the differential equation 

dU int 



71 , 



ih- 



dt 



Utot(t = 0) = 1. 



(57) 

(58) 
(59) 



We want to consider the case where and r/ are on opposite sides of the barrier. In the 
limit when the initial and final states are far away from the barrier, the transition amplitude 
is given by the S-matrix element, which can be expressed in terms of the propagator as 
(Balantekin and Takigawa, 1985) 



Sn f ,nXE) 



— lim 

%h r *- > °° 

r j — ' — oo 



PjP/ 



cxp 



(p/r/ - Pii 



J dTe^ ET l n K{r h n h T ] r h n l ^), 



(60) 



where pj and pf are the classical momenta associated with r*j and ry. In heavy ion fusion 
we are interested in the transition probability in which the internal system emerges in any 
final state. For the Eih partial wave, this is 



T e (E)=Y,\Sn f ,nAE)\' 

which becomes, upon substituting Eqs. ([57]) and (60), 



(61) 



lim 

r f — * — oo 



PiPf 



dTexp 



-ET 



Texp 



ET 

n 



V[r(t)} / X?[f(t)]exp 



-(S(r,T) - S(r,T)) 



p M (r(t),T;r{t),T). 



(62) 



Here we have assumed that the energy dissipated to the internal system is small compared to 
the total energy and taken p f outside the sum over final states. We identified the two-time 
influence functional as 



p M (m,T;r(t),T) = ^ W* ftn ,(r(t);f, 0)W nfini (r(t);T, 0). 

Uf 



(63) 
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Using the completeness of final states, we can simplify this expression to write 



p M (r(t),T;r(t),T) = (n, \uU?®, f WinMt), T)\ m) ■ (64) 

Eq. (|6~4]) shows the utility of the influence functional method when the internal system has 
symmetry properties. If the Hamiltonian in Eq. (^) has a dynamical or spectrum generating 
symmetry, i.e., if it can be written in terms of the Casimir operators and generators of a 
given Lie algebra, then the solution of Eq. ( |5*8"D is an element of the corresponding Lie group 
(Balantekin and Takigawa, 1985). Consequently the two time influence functional of Eq. 
(|64l) is simply a diagonal group matrix element for the lowest-weight state and it can be 
evaluated using standard group-theoretical methods. This is exactly the reason why the 
path integral method is very convenient when the internal structure is represented by an 
algebraic model such as the Interacting Boson Model. 

Two-time influence functionals can be calculated exactly for only a limited number of 
systems. One of these is a harmonic oscillator, linearly coupled to the translational motion. 
In this case the Hamiltonian is 

h 2 d 2 1 
H = + K °( r ) + ( flta + o)^ + «o/(r)(a + a f ), (65) 

-£pQ OT L 

where po is the bare mass of the macroscopic motion and Voif) is the bare potential. The 
(a), lo, m and f(r) are the creation (annihilation) operators of the oscillator quanta, the 
frequency and the mass parameter of the oscillator and the coupling form factor, respectively. 
The quantity «o — [h/2mu] 1 ^ 2 represents the amplitude of the zero-point motion of the 
harmonic oscillator. The two time influence functional pu reflects the effects of coupling to 
the harmonic oscillator and is given by 

~ ~ % ~ Q 2 

p M {r{t),T; r(t),T) = exp[--u;(T - T)\ exp{-^( yi + y 2 + y 3 )} (66) 



with 

2/1 

V2 



T dt f ds f(r(t))f(r(s))e~ iuj{t - s) 

I! JO 

T dt T ds f(r(t))f(m)e Mi - S) 



o Jo 



y 3 = _ e M?-r> [ T dt f(r(t))e^ [ T dt f{f{i))e-^. (67) 
Jo Jo 

An exact calculation of the influence functional is possible in this case because of the sym- 
metry of the Hamiltonian under the Heisenberg-Weyl algebra. 



D. Adiabatic and Sudden Tunneling 

We can discuss the effects of couplings between nuclear structure and translational motion 
in two limiting cases. The first case is the sudden limit in which the energy levels of the 
internal system are degenerate. The second case is the adiabatic limit in which the energy 
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of the first excited state of the system is very large so that the internal system emerges in 
the ground state at the other side of the barrier. 

Several examples are worked out explicitly by Balantekin and Takigawa (1985), to which 
the reader is referred for further details. It can be shown that in the sudden limit, the total 
transmission probability is given by an integral of transmission probabilities for a fixed value 
of the internal coordinate, with the weight of the integration given by the distribution of 
the internal coordinate in its ground state. For example for the linearly-coupled harmonic 
oscillator as the excitation energy gets smaller, u — > 0, the influence functional of Eq. fl66|) 
takes the form 

p M (r(t),f;r(t),T) = exp L^[£ d tf(r(t)) - jTdf/(r(f))] a j . (68) 
Using the integral 

J°° dxe-( ax2+b ^ = y| e b2 ' Aa (69) 

we can rewrite the influence functional as 

p M (r(t),f; r(t),T) = p S (r(t),f; r(t), T) 

l r°° , _i(-s-)2 

dae 2ya o' 



aov 2vr J-co 

x exp dtf(r(t)) - £ dtf(r(t))]j (70) 

where the lower suffix s stands for sudden tunneling. Inserting this expression into Eq. (p2[) 
one obtains 

1 r+oo 

T(E) = -= / dxe- x / 2 T (E, V (r) + xa f(r)), (71) 

V ZTC J-co 

where T (E, Vo(r)+xaf(r)) is the probability of tunneling through the one-dimensional bar- 
rier Vo(r) + xaf(r) at energy E. This expression is known as the zero point motion formula 
and was first derived by Esbensen (Esbensen, 1981). This result can also be derived either 
using the coupled-channels formalism (cf. Section |111 B| ) or using Green's functions (Taki- 
gawa et al, 1992a). Similarly, if the translational motion couples to a very slow rotational 
motion through a coupling Hamiltonian given by 

H int = J^l3P 2 (cos9)f(r) (72) 

V 47T 

then the net tunneling probability is obtained by first calculating the tunneling probability 
for a fixed orientation of the principal axis of the deformed nucleus, and then taking average 
over all orientations 

T{E) = £dcos9T (E,V (r) + J ^-{3P 2 (cos6)f(r)\ . (73) 
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This formula was first derived by Chase et al. (1958) in the study of scattering of rotational 
nuclei in the sudden approximation. Systematics of the fusion cross sections of 16 with a 
series of Sm isotopes, ranging from vibrational to rotational, was first given by Stokstad and 
Gross (1981) using Eq. flTj). 

In actual calculations, both vibrational and rotational excitations are truncated at finite 
excited states. In these cases, the Hermite and the Gauss integrals in Eqs. ([Q]) and (|73|) are 
replaced by the Hermite and the Gauss quadratures, respectively (Nagarajan et al, 1986). 
Zero point motion formulae then have a simple geometric interpretation: In this approxi- 
mation fusion of a deformed nucleus with a finite number (N) of levels can be described by 
sampling N orientations with their respective weights: 

N 

a(E) = J2ui<T{E,V {R) + Xif{r)), (74) 

i=l 

where J2f=i^i — 1- For example, in a two-level system, the orientations 9i = 70.12° and 
6 2 = 30.55° contribute with weight factors uji = 0.652 and uji = 0.348, respectively, as 
illustrated in Figure 0. 

When both target and projectile are heavy and deformed one needs a model describing 
macroscopic potential energy surfaces for arbitrarily oriented, deformed heavy ions. Such 
a model describing completely general configurations of two separated nuclei is given by 
Moller and Iwamoto (1994). 

In the adiabatic limit, i.e., that of slow tunneling, one can introduce (Balantekin and 
Takigawa, 1985) an Vl/u expansion, uj and Q being the frequencies of the internal motion and 
the tunneling barrier, respectively. In this limit the effects of the coupling to internal degrees 
of freedom can be represented in terms of energy-independent potential renormalization 
(Balantekin and Takigawa, 1985; Muller and Takigawa, 1987) : 

V (r) - V ad = V (r) - otlf{r)/hw, (75) 

and a mass renormalization (Takigawa et al., 1994a) 

This limit is very closely related to the situation investigated by Caldeira and Leggett (1983) 
in their study of multidimensional quantum tunneling. In the Caldeira-Leggett formalism 
the term which renormalizes the potential in Eq. ( [75]) is added to the Lagrangian as a 
counter-term. Note that the correction term to the potential is the well-known polarization 
potential and the correction term to the mass is the cranking mass. Hence appropriate 
generalizations of Eqs. ( ff5|) and (|76|) hold not only for a linearly-coupled oscillator, but for 
any system (Takigawa et al, 1994a). Effects of the polarization potential on the subbarrier 
fusion cross section were elucidated by Tanimura et al. (1985). 

Typically adiabatic barrier alone overpredicts the transmission probability. We can 
demonstrate this for the constant coupling case (/i a d = /io) using the exact coupled-channels 
result of Eq. (|5| 

TUE) = f 1 exp(-a S /7(M 2 )T (e - nfU, + £f£) . (77) 
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Inserting the inequality 



TIE — nhu + 




) 



<T\E + 




(78) 



into Eq. ([77]) we get 



T CC (E) < T adiablltic (E) . 



(79) 



We should emphasize that in the adiabatic limit both the potential and the mass renor- 
malization should be considered together. Many authors refer to using only the adiabatic 
potential as the adiabatic limit. While for constant coupling it is true that /i ac j = /io, in 
most cases of interest to nuclear fusion the coupling form factors rapidly change near the 
barrier and the difference between the adiabatic mass and the bare mass can get very large. 
In these cases, using the adiabatic correction to the mass in addition to the adiabatic po- 
tential, even though the former is in the next order in 1/u as compared to the latter, can 
significantly reduce the transmission probability to values below the exact coupled-channels 
result (Takigawa et ai, 1994a). 

Adiabatic and sudden approximations are very useful for obtaining analytical results 
which provide a conceptual framework for understanding the fusion process. For deformed 
nuclei, where the excitation energies are very low, sudden approximation provides a reason- 
ably good description of the data. In case of rotation-vibration coupling, sudden approxi- 
mation can be utilized to reduce the size of the channel coupling. 

If the excitation energies of the internal system are large the sudden approximation tends 
to overestimate the tunneling probability at energies well below the barrier. Indeed Esbensen 
et al. (1983) showed that at below the barrier energies sudden approximation provides an 
upper limit to the tunneling probability. Using the path integral approach Hagino et al. 
(1995b) showed that, for a nearly degenerate system, the finite excitation energy leads to 
a multiplicative dissipation factor which reduces the barrier penetrability estimated in the 
sudden limit. This latter result also provides a good example of the utility of the path integral 
method. As we mentioned earlier for Hamiltonians written in terms of the generators of a 
Lie algebra path integral approach is especially convenient as integrals over paths become 
integrals over the group manifold and thus are amenable to the standard group theoretical 
techniques. An example for this is the treatment of the nuclear structure effects in subbarrier 
fusion with the interacting boson model discussed in Section [TV| Finally coupled channel 
techniques are not practical as the number of channels gets very large, e.g., for tunneling 
problems at finite temperature. On the other hand for finite temperature problems path 
integral techniques can easily be applied (see e.g., Abe and Takigawa (1988) where particle 
decay from a hot compound nucleus is investigated and the thermal fluctuation of the nuclear 
surface is shown to amplify the dynamical effect). 



In the intermediate cases between the sudden and the adiabatic tunneling, the effects of 
the environment is not straightforward to illustrate in simple physical terms. The dynamical 
norm method (Brink, 1985b; Takigawa et ai, 1995b) is a technique introduced to give an 



E. Intermediate Cases - Dynamical Norm Method 
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intuitive understanding of the effects of the environment in the intermediate cases. Defining 
the adiabatic basis by 

[Ho{Z) + H int ((,r{t))]<f> n (r{tU) = h(r(t),0<f> n (r(t),0 = e n (r(t))<f> n (r(t),0- (80) 

In the same spirit as the WKB approximation for a potential model, for the total wave 
function one can take the ansatz 

*(r,fl = *(r,*)(^)" e fc *W. (81) 

Here the parameter e 2 is 1 and —1 in classically allowed and in classically forbidden regions, 
respectively and the action obeys the Hamilton- Jacobi equation 



2 uw\ 2 



+ V ad (r) = E (82) 



2fj, \ dr 

where V a d is the adiabatic potential defined by 

V ad (r) = V (R) + e (r). (83) 

One can then derive an approximate expression where the net tunneling probability is given 
by the product of the tunneling probability through the adiabatic potential barrier and a 
multiplicative factor which represents the non-adiabatic effect of the intrinsic degrees of 
freedom : 

T = T (E,V ad )-Af(r b ) 
A/-(r)= f dt\*{r(TU)\ 2 , (84) 



where r& is the time when the tunneling process is completed. Note that in this approxi- 
mation the transmission probability in Eq. ( [S4]) is calculated with the bare mass, not the 
renormalized mass of Eq. (|76|) . One can easily show that J\f < 1 . This means that the trans- 
mission probability calculated with the adiabatic potential is always greater than the actual 
transmission probability. The deviation of TV from 1 gives the measure of non-adiabaticity 
of the tunneling process. 

One way to determine the dynamical norm factor, A/", is to first expand $(r, £) in the 
basis of the adiabatic states 0„(r(r),£) 

$(r,0=T, a n(T)-<t>n{r{T),0 (85) 

n 

and then determine the expansion coefficients a n by solving the coupled-channels equations 

i " w "' (r) 5:" (I) M^w) 

dh 1 

< 4>n\-^\4> m >= --e n (r(r))a n (T) (86) 
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with i n = e n (r(r)) — e (r(r)). Note that Eq. (|86| ) is correct to describe a classically forbidden 
region. One needs to modify it with e in order to describe a classically allowed region. As 
an example, one can consider the linearly-coupled oscillator with the coupling form factor 
given by /(r) = cr. Assuming that the tunneling path is given by t(r) = i? sin(f2r), Ro 
being the length of the tunneling region, we obtain 

Eq. shows that the adiabaticity of the tunneling process is governed not only by a 

parameter Q/u, but also on the coupling strength and on the length of the tunneling region. 
Takigawa et al. (1995b) applied the dynamical norm method to spontaneous fission of 234 U. 

IV. DESCRIPTION OF NUCLEAR STRUCTURE EFFECTS BY THE 
INTERACTING BOSON MODEL 

An algebraic nuclear structure model significantly simplifies evaluation of the path in- 
tegral. The Interacting Boson Model (IBM) of Arima and Iachello (Iachello and Arima, 
1987) is one such model which has been successfully employed to describe the properties of 
low-lying collective states in medium-heavy nuclei. In this section, attempts to use the IBM 
in describing nuclear structure effects on fusion are reviewed. The path integral formulation 
of this problem, as sketched in the next section, requires analytic solutions for the nuclear 
wave functions. In the first attempt to use the IBM to describe nuclear structure effects 
in subbarrier fusion, the SU(3) limit of the IBM was employed (Balantekin et al, 1991). 
However, the SU(3) limit corresponds to a rigid nucleus with a particular quadrupole defor- 
mation and no hexadecapole deformation, a situation which is not realized in most deformed 
nuclei. Thus analytic solutions away from the limiting symmetries of the IBM are needed 
for realistic calculations of subbarrier fusion cross sections. 

In a parallel development, a 1/N expansion was investigated (Kuyucak and Morrison, 
1988; Kuyucak and Morrison, 1990; Kuyucak and Morrison, 1991) for the IBM which pro- 
vided analytic solutions for a general Hamiltonian with arbitrary kinds of bosons. This 
technique proved useful in a variety of nuclear structure problems where direct numerical 
calculations are prohibitively difficult. Later it was applied to medium energy proton scat- 
tering from collective nuclei (Kuyucak and Morrison, 1989) in the Glauber approximation, 
generalizing the earlier work done using the SU(3) limit (Ginocchio et al., 1986). Using the 
1/N expansion technique in the path integral formulation of the fusion problem (Balantekin 
et al, 1992) makes it possible to go away from the three symmetry limits of the IBM. In 
particular arbitrary quadrupole and hexadecapole couplings can be introduced. 



A. Linear Coupling 



We take H int in Eq. (|33| ) to be of the form of the most general one-body transition 
operator for the IBM, 



rat 



b% 



(k) 



r<*>(r), 



(88) 
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where the boson operators are denoted by bi and &j. The k sum runs over k — 2, 4, . . . 2£ max . 
Odd values of are excluded as a consequence of the reflection symmetry of the nuclear 
shape, and the k = term is already included in the bare potential Vo(r). The form 
factors akji(r) represent the spatial dependence of the coupling between the intrinsic and 
translational motions. The interaction term given in Eq. (|88"D is an element of the SU(6) 
algebra for the original form of the Interaction Boson Model with s and d bosons and is 
an element of the SU(15) algebra when g bosons are included as well (Iachello and Arima, 
1987). 

To simplify the calculation of the influence functional, we use the no-Coriolis approxima- 
tion. We first perform a rotation at each instant to a frame in which the z-axis points along 
the direction of relative motion. Neglecting the resulting centrifugal and Coriolis terms in 
this rotating frame is equivalent to ignoring the angular dependence of the original Hamil- 
tonian. In this approximation, the coupling form factors become independent of I and only 
m = magnetic sub-states of the target are excited (Takigawa and Ikeda, 1986; Tanimura, 
1987; Esbensen et al, 1987). For heavy systems the neglected centrifugal and Coriolis forces 
are small. We take the scattering to be in the x-y plane. Then making a rotation through 
the Euler angles b = (0, 7r/2,0), we can write the full Hamiltonian as the rotation of a 
simpler Hamiltonian depending only on the magnitude of r 

H = R(b)H(°\r)rf(b). (89) 

Since in Eq. (|3"B|) H Q and H t = + Vq(t) are rotationally invariant, H int is the only term 



whose form is affected by the transformation. Hence we introduce the rotated interaction 
Hamiltonian (r), given by 

H int = R(h)H^(r)rf(b), (90) 
H^(r) = J2Mr)b] m bim, (91) 



jlm 



9 j£m { 



2k 4- 1 

E V—j^-ti™* ~ m\k0)a kji (r). (92) 



k 



If we assume now that the form factors akje(r) are all proportional to the same function of 
r then the Hamiltonian commutes with itself at different times and hence we can write 
the two-time influence functional as 

pT \ / 7' rT 



Pm = ( rii 



11; 



(93) 



in the degenerate spectrum limit. 

Since the exponents of the two operators in the influence functional commute, pu be- 
comes the matrix element of an SU(6) transformation between SU(6) basis states, in other 
words it is a representation matrix element for this group and can easily be calculated us- 
ing standard techniques. The two-time influence functional for the sd-version of the IBM 
was calculated by Balantekin et al. (1992) and, for the particular case of SU(3) limit, by 
Balantekin et al. (1991). 
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B. Higher-Order Couplings 



Up to this point, we have utilized only a first-order coupling between nuclear states 
and translational motion. Alternatively, one can include the effects of coupling to all orders. 
This can be achieved by exploiting the symmetry properties of the resolvent operator directly 
without utilizing its path integral representation. Such a Green's function approach has also 
been used to study quantum tunneling in a heat bath (Takigawa et al, 1992a). 

To include the effects of couplings to all orders, the interaction Hamiltonian in Eq. ( |33"D 
is written as 

H mt (r, + V {r) = V Coul (r, f) + V nuc (r, 0, (94) 
where the Coulomb part is 

Mr,0~(l + ^O) (r>R 1 ), 
r 5 r A 

= ^(l + jjgd) (r<ft). (95) 
The nuclear part is taken to have Woods-Saxon form, 

K,(r, = -V (l + exp ( T ~ "» " Rl6{i ■ ) ) " . (96) 



In Eqs. ( j9q) and fl96|), i?o is the sum of the target and projectile radii and R\ is the 
mean radius of the deformed target. O is a general coupling operator between the internal 
coordinates and the relative motion 

6 = ]T^T( fc )(0-r«(r). (97) 

k 

The coefficients V}. represent the strengths of the various multipole transitions in the target 
nucleus. In the standard IBM with s and d bosons, the only possible transition operators 
have k — 0, 2, 4, . . . (odd values being excluded as a consequence of the reflection symmetry 
of the nuclear shape). The monopole contribution is already included in the Woods-Saxon 
parameterization and so is not needed. The quadrupole and hexadecapole operators are 
given by 

= [ s t d ~ +d t s p) +x [Jdf\ (98) 
= [rft^W. (99) 

We adopt the "consistent-Q" formalism of Casten and Warner (Casten and Warner, 1988), 
in which x i n Eq. (|98| ) is taken to be the same as in H IBM (fitted to reproduce the energy 
level scheme and the electromagnetic transition rates of the target nucleus) and is thus not 
a free parameter. 

In the previous section, we used the usual approximation in which the nuclear potential 
of Eq. (PB|) is expanded in powers of the coupling, keeping only the linear term (cf. Eq. 
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(|38|)). In order to calculate the fusion cross section to all orders we consider the resolvent 
operator for the system 

G+(E) = ^. (100) 

E+ - H t - H IBM (Z) - H int (r,0) 

The basic idea is to identify the unitary transformation which diagonalizes the operator O 

O d = UOU ] (101) 

in order to calculate its eigenvalues and eigenfunctions 

O d \n) = Cn\n). (102) 
Assuming the completeness of these eigenfunctions 

^2\n)(n\ = l (103) 

n 

one can write the matrix element of the resolvent as 

( t f ,r f \G + {E)\t i ,r i ) 



{i f ,r f \U ] U [E+ -H t - H int (r,6)\ Vwi&n) 
(Z f ,r f \rt [E + -H t - H- mt (r, O d )] 1 £ \n) (n|W|&, r<> 

n 

E^/l^lnXnlWlftXr/IG+Iri) (104) 



where 



G ^> = E+-H t -H M <105) 

To derive this result we ignore the excitation energies in the target nucleus. This corresponds 
to setting the term H IBM to be zero in Eqs. (|100|) and ( |105|) . In this case, the G+(E) given 



in Eq. ( |104j ) is the resolvent operator for one-dimensional motion in the potential H- mt (r, £ n ), 
the fusion cross section of which can easily be calculated within the standard WKB approx- 
imation. The total cross section can be calculated by multiplying these eigen-channel cross 
sections by the weight factors indicated in Eq. ( |104|) . The calculation of the matrix ele- 
ment (n\U\^i) within the IBM is straightforward (Balantekin et ai, 1993; Balantekin et al, 
1994a; Balantekin et ai, 1994b). It is also possible to generalize the previous formalism to 
include arbitrary kinds of bosons in the target nucleus and investigate whether g bosons 
have any discernible effects on subbarrier fusion reactions. One finds that (Balantekin et 
ai, 1994b) except for slight differences in the barrier distributions (which can be made even 
smaller by fine tuning the coupling strengths), there are no visible differences between the 
sd and sdg model predictions. The similarity of the results implies that subbarrier fusion 
probes the overall coupling strength in nuclei, but otherwise is not sensitive to the details 
of the nuclear wave functions. In this sense subbarrier fusion, which is an inclusive process, 
is in the same category as other static quantities (energy levels, electromagnetic transition 
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rates), and does not seem to constitute a dynamic probe of nuclei, in contrast to exclusive 
processes such as proton scattering. 

Using this formalism, a systematic study of subbarrier fusion of 16 with rare earth nuclei 
became possible. Fusion cross sections for the reactions 16 + 144 > 148 > 154 Sm and 16 + 186 W 
were measured by the Australian National University group (Wei et ai, 1991; Leigh et al, 
1993; Lemmon et ai, 1993; Morton et al, 1994; Leigh et ai, 1995). The angular momentum 
distributions for 16 + 154 Sm was measured by Bierman et al. (1993), and for 16 + 152 Sm by 
Wuosmaa et al. (1991). Those for 16 + 144,148 Sm were deduced by Baba (1993). Balantekin 
et al. (1994a) fit the existing data on vibrational and rotational nuclei with a consistent set 
of parameters which are then used to predict the cross-section and < i >-distributions. 
Figure |8| compares those data with the cross section calculations of Balantekin et al. (1994a) 
and the angular momentum distribution calculations of Balantekin et al. (1994c). 

Finally one should mention that the effects of anharmonicities both in nuclear spectra 
(Casten and Zamfir, 1996) and in the vibrational coupling in subbarrier fusion (Hinde et ai, 
1994; Morton, 1995) recently attracted some attention. Arima and Iachello (1976) pointed 
out that the U(5) symmetry limit of the IBM should exhibit in the spectra anharmonicities 
similar to the geometric anharmonic vibrator model of Brink et al. (1965). This assertion 
was later explicitly confirmed (Aprahamian et ai, 1987). This feature of the IBM makes it 
possible to discuss the effects of anharmonicities in the vibrational coupling in subbarrier 
fusion using the Z7(5) limit (Hagino et ai, 1997b). 

V. COMPARISON OF CURRENT THEORY WITH DATA 
A. Status of Coupled-Channels Calculations 

DiGregorio and Stokstad (1991) presented a systematic analysis of fusion cross sections 
and average angular momenta for fourteen different systems using a barrier penetration 
model that includes coupling to inelastic channels. They concluded that model predictions 
explain data well for light and asymmetric systems whereas large discrepancies exist for 
large symmetric systems. Indeed for light and asymmetric systems the basic premise of the 
coupled-channels calculations is justified: For these systems the repulsive Coulomb potential 
is relatively weak and the tail of the attractive nuclear potential has sufficient strength to 
"turn it around" to form the potential barrier. The barrier is thus formed at a rather 
large nuclear separation, long before two nuclear surfaces start touching. Consequently, as 
the system penetrates the barrier individual nuclei preserve their character and one can 
talk about coupling of the states in the target nuclei to the quantum tunneling process. 
The inversion procedure of Balantekin et al. (1983) demonstrated that even when there are 
isotopic differences in the fusion cross section (Wu et al, 1985), one can still describe the 
quantum tunneling with a one-dimensional effective potential for very light systems. On the 
other hand, for heavier and more symmetric systems other effects not explicitly included in 
the coupled-channels calculations, such as neck formation may play an important role. 

Many experimental groups use simplified coupled-channels codes such as CCFUS (Dasso 
and Landowne, 1987b), CCDEF (Fernandez-Niello et al, 1989), CCMOD (Dasgupta et al, 
1992), or the IBM-based models (Balantekin et al., 1993; Balantekin et al., 1994a; Balantekin 
et al., 1994b). As we mentioned earlier it is worthwhile to keep in mind that although these 
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codes are quite adequate for qualitative comparisons, they may be making a number of 
assumptions, such as ignoring the radial dependence of the coupling form factor, excitation 
energies and/or higher-order couplings. Before one makes a quantitative statement, it may 
be sagacious to check what the approximations are and if the observed discrepancies with 
the data is a result of these simplifications. In Section [IV B, we explicitly demonstrated the 
effect of the higher order couplings. 

The crucial ingredient of the coupled-channels calculations is to identify relevant degrees 
of freedom and to model the appropriate Hamiltonian. Even the choice of the optical 
potential should be scrutinized. For example, a recent survey (Brandan and Satchler, 1997) 
of the knowledge of the optical potential between even much simpler systems, such as two 
light ions, indicate that many anomalies needed to be resolved before a good theoretical 
understanding of the elastic scattering data can be achieved. For heavier systems, from fits 
to elastic scattering data at energies near the barrier, the optical potential was shown to have 
a strong energy dependence, known as "threshold anomaly" (Satchler, 1991). Nagarajan et 
al. (1985) pointed out that the dispersion relation between the real and imaginary parts 
of the optical potential should be used in regions where the absorption varies rapidly with 
energy such as near and below the barrier. The source of the energy dependence can be 
either channel coupling or the non-locality of the exchange contribution (Galetti and Candido 
Ribeiro, 1994). Galetti and Candido Ribeiro (1995) compared non-local effects and coupled 
channels calculations in simple models of nuclear fusion. 

The very first coupled-channels calculations for heavy-ion fusion assumed a linear cou- 
pling to quadrupole or octupole surface vibrations and quadrupole deformations. As more 
precise data became available, the significance of the hexadecapole deformations (Rhodes- 
Brown and Oberacker, 1983), neutron transfer (Broglia et al., 1983a), coupling of multi- 
phonon states (Takigawa and Ikeda, 1986; Kruppa et al, 1993), and higher-order couplings 
(Balantekin et al, 1993) emerged. In the rest of this section we discuss representative data 
illustrating these effects. 



B. Nucleon Transfer 

Another interesting question is the effect of the nucleon transfer on subbarrier fusion (cf. 
Section |111B| ). In particular the role of transfer channels with positive Q- values has been 
emphasized (Broglia et al., 1983a; Broglia et al., 1983b; Esbensen and Landowne, 1989). 
Algebraic models, such as that described in Section [TVL at present do not include the effects 
of nucleon transfer. It is now experimentally possible to observe up to six-nucleon transfer 
at subbarrier energies (Jiang et al, 1994). Hence in the near future systematic studies 
including nucleon transfer reactions, fusion, and elastic scattering may be possible. 

The effect of nucleon transfer on fusion can be illustrated for example by considering 
fusion reactions between different Ni isotopes. Indeed these were the pioneering experiments 
of Beckerman et al. (1980) where the enhancement of subbarrier fusion cross sections was 
first observed. These cross sections were later measured by Schicker et al. (1988) and more 
recently by Ackermann et al. (1996). The cross sections measured by Ackermann et al. (1996) 
for the 58 Ni+ 64 Ni and 64 Ni+ 64 Ni systems are displayed in Figure |9|, where the energies are 
normalized to the height of the s-wave potential barrier. One sees a discernible enhancement 
for the 58 Ni+ 64 Ni system over the 64 Ni+ 64 Ni system. In the symmetric system 64 Ni+ 64 Ni 
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there are no transfer channels with positive Q-values, and only those channels describing 
inelastic excitations need to be included. On the other hand, in the 58 Ni+ 64 Ni system 
there is an additional coupling of the transfer channel 64 Ni( 58 Ni, 60 Ni) 62 Ni with a Q-value 
of Q = +3.9 MeV. These additional channels increase the 58 Ni+ 64 Ni cross section (cf. the 
discussion in Section |111 B| ) . 

Signatures of positive Q-value transfer reactions can also be identified in fusion barrier 
distributions. By comparing barrier distributions for 16 0+ 144 Sm and 17 0+ 144 Sm reactions 
Morton et al. (1994) showed that the effect of the neutron-stripping channel in the second 
reaction is evident in the barrier distribution. 

In fusion reactions of identical nuclei there are a number of interesting effects magnified 
by the existence of only even partial waves. For example, the fusion cross sections have an 
oscillatory structure as a function of energy (Poffe et al., 1983). Furthermore, elastic transfer 
plays an important role in such collisions (von Oertzen and Noremberg, 1973). Christley et 
al. (1995) showed that even in cases where no oscillatory structure is visible in cross section, 
there still is a signature of the elastic transfer in the barrier distributions. 



C. Probing Asymmetry Effects 

One way to probe the effects of asymmetry of the system other than nucleon transfer is 
to measure fusion cross sections and average angular momenta for different systems leading 
to the same compound nucleus. Such a measurement was recently performed (Ackermann 
et al., 1996) for the systems 28 Si+ 100 Mo and 64 Ni+ 64 Ni leading to the compound nucleus 
128 Ba. These measurements complement a previous measurement for the system 16 0+ 112 Cd 
(Ackermann et al., 1994). They find that both fusion cross sections as well as average 
angular momenta can be explained by coupled channels calculations. For the 28 Si+ 100 Mo 
system including lowest 2 + and 3~ states of both the target and the projectile in the coupled- 
channels calculation improves the agreement between theory and data, as compared to the 
no-coupling limit, but is not sufficient to reproduce the data. One needs to include an 
additional channel with the Q-value of the two-neutron pick-up reaction to bring the data 
and theory into agreement. The data for the symmetric system 64 Ni+ 64 Ni, where no transfer 
channels with positive Q-values are present, are already well reproduced only with coupling 
to the inelastic channels. Hence the data of Ackermann et al. (1996) also provide evidence 
for the influence of the two-nucleon transfer channels with positive Q-values on the fusion 
probabilities. 

Studies of transfer channel coupling and entrance channel effects for the near and subbar- 
rier fusion was also given by Prasad et el. (1996) for the systems 46 Ti + 64 Ni, 50 Ti + 60 Ni, 19 F 
+ 93 Nb, and by Charlop et al. (1995) for the systems 28 Si+ 142 Ce, 32 S+ 138 Ba, and 48 Ti+ 122 Sn. 
These authors also report no significant entrance channel effects except that the positive Q 
value for two-neutron pickup shows up as an additional enhancement in the 46 Ti + 64 Ni 
system. 
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D. Signatures of Nuclear Vibrations 

One relatively unexplored aspect of subbarrier fusion is searching for signatures of nuclear 
vibrations. Leigh et al. (1995) confirmed the effects of vibrational coupling in the 16 0+ 144 Sm 
system. To search for signatures of nuclear vibrations Stefanini et al. (1995c) measured the 
fusion cross section for 32 . 36 s+ 110 Pd systems. 110 Pd is a vibrational nucleus whose two- 
quadrupole-phonon triplet is well known. Simplified coupled-channels calculations for these 
systems are performed by Stefanini et al. (1995c) and by DeWeerd (1996). Stefanini et al. 
(1995c), using the method of Kruppa et al. (1993), assumed a constant coupling explicitly 
including the finite Q value of the coupled channels. In these experiments it was also 
observed that the cross section for the 32 S+ 110 Pd system is greatly enhanced because of the 
two-neutron transfer channel. 

Since 110 Pd lies between U(5) and SO(Q) symmetry limits of the IBM it cannot be 
described analytically. DeWeerd (1996) first numerically calculated quadrupole matrix ele- 
ments between different states using the PHINT code (Scholten, 1991) and then numerically 
obtained eigenvalues and the associated weights (cf. Eq. ( |104j )). His result for the barrier 
distribution of the 36 S+ 110 Pd system is displayed in Figure |10] along with the data. While 
one-phonon space clearly fails in describing the barrier distribution, the agreement with data 
successively improves as one includes more phonons in the calculation. This drastic change 
in the barrier distribution for different numbers of phonons was also noted by Stefanini et 
al. (1995c) for this system and by Stefanini et al. (1995b) for the 58 Ni+ 60 Ni system. 



E. Effect of Non-linear couplings 

An important component in the theoretical description of the subbarrier fusion data is 
the effect of non-linear couplings. These effects were discussed by Balantekin et al. (1993) 
using the IBM in the limit of zero excitation energy (cf. Section [I V B|) . For nuclear surface 



vibrations the excitation energies cannot be neglected in most cases, and one has to solve 
full coupled-channels equations. These calculations can be very involved, and consequently 
they were carried out by very few groups. 

Esbensen and Landowne (1987) expanded the coupling potential (cf. Eq. (0)) up to 
second order with respect to the deformation parameter obtaining a good agreement be- 
tween their calculations and data for fusion cross sections between different nickel isotopes. 
The quadratic coupling approximation was shown to describe well fusion cross sections and 
angular momentum distributions for the 58 > 64 Ni+ 92 ' 100 Mo reactions (Rehm et al, 1993). 
Coupled-channels calculations including coupling to all orders and the finite excitation en- 
ergy of nuclear surface vibrations were performed for the 58 Ni+ 60 Ni reaction (Stefanini et 
al., 1995b). 

In Figure [TT] the calculation of Hagino et al. (1997a) for the system 64 Ni+ 96 Zr is compared 
with the data of Stefanini et al. (1992). Here the results of coupled-channels calculations to 
all orders (solid lines) agree with the data very well as opposed to the no-coupling (the dotted 
lines), linear-coupling (the dot-dashed lines), and quadratic coupling (the dashed lines) cases. 
The upper panel compare theory and calculations for the fusion cross section and the lower 
panel for the average angular momenta. In this calculation up to double phonon states 
are included in the coupled channels. An important feature of this calculation is that a 
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truncation of the coupling within the double phonon space even at the quadratic level is not 
sufficient to describe the data, one needs to include couplings to all order. 



F. Angular Momentum Distributions 

DiGregorio and Stokstad (1991) also compared average angular momenta obtained using 
different experimental techniques with theoretical predictions. They found that there is good 
agreement between theory and data obtained from isomer ratio and gamma-ray multiplicity 
measurements with the exception of more symmetric systems, but not for the fission fragment 
measurements. From the fission fragment angular anisotropy measurements one obtains not 
< £ >, but < £ 2 >. Especially if the distribution is pushed to higher £ values as a result of 
coupling to other channels < £ 2 > 1//2 may significantly differ from < £ >, which may explain 
some of the reported discrepancies. An excellent review of the efforts to measure angular 
momentum distributions in fusion reactions was given by Vandenbosch (1992). 

There have been many attempts to extract average angular momenta directly from the 
fusion excitation functions (Balantekin and Reimer, 1986; Reisdorf et ai, 1985; Rowley 
et al, 1993; Balantekin et ai, 1996). Moments of angular momenta are related to the 
moments of fusion cross sections (Balantekin et al, 1996). These relations can help assess 
the consistency of accurate fusion cross section measurements with measurements extracting 
angular momentum distributions using different methods. 



G. Probing Shape Phase Transitions with Fusion 

Until recently little attention was paid to subbarrier fusion on gamma unstable targets. 
The Os and Pt region is interesting to study since these nuclei go through a shape transition 
from prolate to oblate as one increases the number of protons from 76 to 78. 192 0s has 
a positive (prolate) quadrupole deformation parameter and a negative hexadecapole defor- 
mation parameter. 194 Pt has a quadrupole deformation parameter similar in magnitude to 
those of 192 Os, but with a negative sign (oblate) and a hexadecapole deformation parameter 
comparable to that of 192 Os in sign and magnitude. The isotopes are similar in all respects 
other than the fa sign. The effect of this shape phase transition on the barrier distributions 
would be noticeable by the skewness toward higher energies for prolate nuclei and toward 
lower energies for oblate nuclei (Balantekin et al, 1994a). The barrier distributions calcu- 
lated using the IBM based model of Section |IV] by Balantekin et al. (1994a) is shown in 



Figure 12 



To understand this effect of the shape phase transition in the barrier distributions, the 

fusion cross sections for transitional nuclei Pt and Os were recently measured using an 16 O 

beam by the Legnaro group (Stefanini et ai, 1995a) and a 40 Ca beam by the Seattle group 

(Bierman et al, 1996a). (Since 40 Ca is a heavier projectile one expects this effect to be 

enhanced). The total fusion cross sections calculated by Balantekin et al. (1994a) for the 
i6 0+ i94p t system 

agree very well with the evaporation residues measured by Stefanini et 
al. (1995a) once the fission cross sections estimated by the statistical model calculations are 
subtracted. 
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The 40 Ca+ 194 Pt and 40 Ca+ 192 Os cross sections measured and the associated barrier 
distributions extracted by Bierman et al. (1996a) are shown in Figure 13, where the results 
are also compared with the CCDEF calculations. The calculations take the excitations of 
the target nucleus into account within the rotational model including both quadrupole and 
hexadecapole deformations. They also take into account the excitation of the projectile to 
the 3~ state at 3.7 MeV and the two neutron transfer reactions from the target nucleus to the 
ground state of 42 Ca. The constant coupling approximations have been used for vibrational 
excitation of the projectile and the transfer reactions, while the radial dependence of the form 
factor of the collective model was used for rotational excitations. The predicted skewness 
of the barrier distribution toward higher energies for the prolate nucleus 194 0s and toward 
lower energies for the oblate nucleus 194 Pt is observed. The second peak for both systems 
is due to the excitation of the octupole state in the projectile. The barrier distributions 
for both systems also exhibit a tail at the lower energies, which is not reproduced by the 
CCDEF calculations. This problem was recently associated with the contribution of the 2n 
transfer reactions to the first excited 2 + state of 42 Ca (Bierman et al., 1996b). 

One should point out that fusion barrier distribution extracted from the 16 0+ 186 W data 
of Lemmon et al. (1993) has the shape expected for a target nucleus with a negative hexade- 
capole deformation. There are pronounced differences between this distribution and that of 
16 0+ 154 Sm (Wei et al., 1991), which are just those expected from a change in sign of the 
hexadecapole deformation. These data hence demonstrate s strong sensitivity of fusion to 
the hexadecapole deformation. 



H. Difficulties in Extracting Barrier Distributions 

Since barrier distributions include the second energy derivative of the cross section, very 
accurate measurements of excitation functions at closely spaced energies are required. Even 
with very high precision data, smooth barrier distributions can only be obtained under 
certain model dependent assumptions such as fine-tuning the energy spacing for calculating 
second derivatives (Izumoto et el, 1995; DeWeerd, 1996). Krappe and Rossner (1995) 
suggested to use integrals, rather than the derivatives of the fusion data to improve model 
independence in the analyses. Unfortunately, moments of the cross sections are even more 
featureless than the cross sections themselves (Balantekin et al., 1996). In contrast, one of 
the main advantages of using barrier distributions is that they bring out important features 
in the data. However, moments of the cross section could be useful as they are related 
to the moments of angular momenta under certain assumptions (Balantekin and Reimer, 
1986; Dasso et al., 1986; Esbensen and Landowne, 1987; Balantekin et al., 1996). 

Even with a very high precision in the fusion cross section it is rather difficult to extract 
fusion barrier distributions at higher energies as the cross section changes very slowly and 
the errors on the barrier distribution grow with energy. To illustrate the reason for this 
behavior consider a set of fusion data measured at a fixed energy spacing AE. The second 
derivative may be approximated by the point-difference formula (Rowley, 1992) 

-lW) n = (KW- ■ (106) 

If the statistical errors on the cross section are a fixed fraction of their measured values 
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Hence the error on the second derivative increases as the cross section increases, whereas 
the second derivative itself gets smaller at higher energies. Furthermore AE must be small 
enough to resolve any interesting structure, which also contributes to the large errors at high 
energies. 

It has been suggested that information about the barrier distributions may be contained 
in the quasi-elastic scattering excitation functions at backward angles (Kruppa et al., 1993; 
Andres et al, 1988). Timmers et al. (1995) recently developed a method to extract a 
representation of the fusion barrier distribution from quasi-elastic excitation functions. They 
found that although this representation of the quasi-elastic scattering data indeed shows the 
general features of the fusion barrier distributions, its sensitivity is reduced at high energies. 
More recently Rowley et al. (1996) showed that the effects of strong coupling are present 
in the barrier distributions from the elastic scattering, but are smoothed out since different 
eigen-barriers have phase differences. Furthermore the effects are also smoothed by weak 
couplings, which appear in first order in the elastic scattering cross section, but only in 
second order in fusion cross section. It would be important to treat the phase problem 
properly in order to obtain information on barrier distributions from the elastic scattering 
data. 



Heavy-ion fusion reactions induced by a halo nucleus or by an unstable neutron rich 
nucleus are very intriguing current subjects of nuclear physics (Ishihara et el., 1993). Several 
groups have performed experiments to examine whether the fusion cross section in such cases 
is significantly different from that in heavy ion collisions induced by the corresponding stable 
isotopes. Yoshida et al. (1995) studied the fusion reactions of n > 10 > 9 Be with 209 Bi at energies 
near the Coulomb barrier. They observed no significant difference in the excitation function 
for the collision of n Be from that of 10 Be. On the other hand, Fukou-Youmbi et al. (1994) 
have reported that the induced fission cross section near the Coulomb barrier is much larger 
in the n Be + 238 U reaction than that in the 9 Be + 238 U reaction. This is a very interesting 
result, though it is not clear yet whether the fission took place via a compound nucleus 
formation, suggesting an enhanced fusion cross section in the case of unstable isotope. 

Takigawa and Sagawa (1991) and Takigawa et al. (1992b) suggested that the fusion 
cross section will be significantly enhanced if one uses a halo nucleus as the projectile. In 
deriving this conclusion, the existence of a stable soft dipole oscillation of the core nucleons 
against the halo neutrons was assumed. Though it is not yet completely settled down, the 
existence of a physical soft dipole oscillation in light halo nuclei is unlikely (Sagawa et al, 
1995). This might explain why the fusion data of Yoshida et al. (1995) do not show any 
characteristically different behavior for the case of n Be projectile. As has been shown by 
Takigawa and Sagawa (1991), the neutron halo itself can enhance the fusion cross section by 
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statically lowering the fusion barrier. This effect alone is, however, not so drastic. Moreover, 
n Be has less pronounced halo property than n Li. 

Using time-dependent Hartree-Fock theory Kim et al. (1994) showed that nucleon trans- 
fer is enhanced for fusion reactions between a stable and an unstable nucleus with neutron 
halo. 

Some debate exists concerning the effects of break up of the halo nuclei. Takigawa et 
al. (1993) have shown that although the large enhancement of the fusion cross section is 
moderated by the break-up, the halo nucleus n Li still leads to a larger fusion cross section 
than the other Li isotopes. One should point out that this conclusion relies on the assumption 
that there exists a soft dipole resonance in n Li. On the other hand, Canto et al. (1995) 
argued that the fusion cross section of a halo nucleus will be hindered by the break-up effect. 
To the contrary, Dasso and Vitturi (1994) contended that the break-up channel enhance the 
fusion probability. In order to reach a definite conclusion, one needs to know more about the 
radial dependence of the break-up form factor and has to treat both the real and the virtual 
break-up precesses in a consistent way including the associated potential renormalization. 

One interesting subject to be explored is the effects of the bond formation due to halo 
neutrons on the fusion cross section. Bertulani and Balantekin (1993) studied this problem 
by using the fusion between 11 Li, consisting of 9 Li core and halo di-neutrons, and 9 Li as an 
example . As n Li and 9 Li approach each other, there is a particular separation distance 
at which both the Coulomb and nuclear potentials between the two cores are small, but 
the two neutrons in the halo can be shared by both cores. One can then investigate the 
effect of this molecular bonding on the fusion of n Li and 9 Li. These preliminary calculations 
indicated (Bertulani and Balantekin, 1993) a very significant enhancement of the fusion cross 
section due to molecular bonding. Takigawa et al. (1995a) showed that this cross section 
is somewhat reduced when different initial conditions are used. The calculations of both 
groups were done in the adiabatic approximation, which tends to overestimate the effect. 
Existence of molecular bonding, and its effect on the fusion process remains to be an open 
question. 

In this connection, we wish to remark that a polarization of the wave functions of the 
valence neutrons, i.e., the admixture of higher orbits, is essential in order for the bond effect 
to be significant. Imanishi and von Oertzen (1995) studied the bond effect in the fusion 
between n Be and 10 Be and showed that a significant polarization of the valence nucleons 
starts to take place even where the core nuclei are still far apart if the binding energy 
of the valence nucleons is small, and that consequently there exists a large bond effect in 
this system. Notice that n Be has two bound states lpf/ 2 an d 2s^, 2 and one low lying 
resonance state ldt, 2 . The hybridization of these opposite parity configurations causes a 
large polarization. 

It is possible that the response of the projectile in heavy ion collisions induced by a 
neutron rich unstable nucleus can be formulated in terms of the coupling to a resonance 
state. Several models exist to calculate the effect of the width of the resonant state on 
quantum tunneling (Balantekin and Takigawa, 1985; Hussein et al, 1995). 
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J. Fusion of Very Massive Systems and Superheavy Nuclei 



As the compound nuclei formed in the fusion get heavier, fission becomes an increasingly 
important de-excitation channel. For such reactions, evaporation residues and fission frag- 
ments must be added to obtain the total fusion cross section. As systems get more massive 
(ZpZt > 1000), fusion starts competing with other reaction channels representing a signifi- 
cant exchange of energy, charge, mass, and angular momenta. Understanding the dynamics 
of fusion and competing reactions for very massive systems is essential, among other things, 
to assess the conditions for the formation of super-heavy elements. These topics are covered 
in a recent review by Reisdorf (1994). 

A significant difference between the fusion of massive nuclei and the fusion of medium 
weight nuclei which we discussed so far is that the fusion cross section for very massive 
systems is not enhanced, but rather hindered. This situation is encountered when the 
product of the atomic numbers of the projectile and target exceeds about 1800 (Reisdorf, 
1994). The incident energy has to be considerably higher than the fusion barrier expected 
from the Bass potential (Bass, 1974) which was determined to fit the fusion data above 
the barrier for medium weight systems with Z p Z t = 64 — 850 in a potential model. This 
excess energy is called the extra push energy. Bjornholm and Swiatecki (1982) attributed 
this hindrance phenomena to the fact that the fission barrier for massive systems locates 
well inside the potential barrier in the entrance channel, and introduced the concepts of 
the extra push and the extra extra push. The former is the energy needed to overcome the 
conditional saddle, i.e., the saddle under the constraint of mass asymmetry, while the latter 
is the energy needed to carry the system inside the unconditional saddle for fission. Though 
there have been quite a number of experimental as well as theoretical studies of the extra 
push energy, its origin and the dependence on various parameters of the system, such as the 
effective fissility parameter in the entrance channel, are not fully understood. One should 
note that the decrease of the fusion cross section with decreasing bombarding energy in 
massive systems, where there exists an extra push, is also much slower than that expected 
in the potential model where there exists an extra push. This indicates the existence of a 
kind of enhancement mechanism of the fusion cross section in massive systems as well once 
the hindrance effect associated with the extra push problem is isolated. 

An interesting problem concerning the fusion of massive nuclei is the synthesis of super 
heavy elements. A significant advance occurred when Hofmann et el. (1996) at GSI working 
with the SHIP velocity filter has succeeded in producing the super heavy element Z=112 by 
the so called cold fusion method using 70 Zn + 208 Pb reaction. This is the heaviest element 
recorded to date and it is only two atomic number away from the predicted magic number 
Z=114. Though the successful synthesis of element 112 after the synthesis of elements 110 
and 111 in 1994 seems to indicate that a similar experimental strategy can be used to go 
further to the realm of the heaviest elements, a problem is that the cross section is very small, 
i.e., of the order of 1 pb. Actually, only two events were identified for Z — 112. It would 
certainly be very interesting to look for alternative ways to synthesize super heavy elements. 
An interesting question in this connection is to examine whether there are advantages of 
using neutron rich unstable nuclei. A preliminary study in this direction has been undertaken 
by Takigawa et al. (1992b) and Takigawa and Shinozuka (1992c). These authors discussed 
the advantages such as the larger survival probability of the compound system, lowering of 
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the fusion barrier, and the possible lowering of the extra-push energies, and disadvantages 
such as the low beam intensity in reactions induced by neutron rich unstable nuclei. In 
passing, we wish to mention that Nomura et al. (1995) is trying to use (HI, axn) reactions 
to experimentally synthesize super heavy elements using the cooling mechanism by a particle 
emission, and that Aritomo et al. (1997) are introducing a diffusion model to theoretically 
discuss the mechanism of the synthesis of super heavy elements, though both of them treat 
a thermal process rather than a quantum tunneling process. 

VI. OPEN PROBLEMS AND OUTLOOK 

Although there still are many unsettled issues even in fusion reactions with stable nuclei, 
remarkable progress has been made in the last fifteen years. New, conceptually alluring 
ways of analyzing data, such as studying barrier distributions; new approaches to channel 
coupling, such as the path integral and Green function formalisms; and alternative meth- 
ods to describe nuclear structure effects, such as those using the Interacting Boson Model, 
were introduced. The roles of nucleon transfer, higher-order couplings, and shape-phase 
transitions were elucidated. We can now understand the data for clean (i.e., asymmetric) 
systems in terms of inelastic excitations and nucleon transfer. Acquisition of high precision, 
complementary data for fusion, transfer reactions, and elastic scattering at below and near 
the barrier should be encouraged as theoretical tools are available to analyze them. 

On the other hand, fusion cross sections for very heavy symmetric systems cannot be 
reproduced by the present models. In such systems inclusion of higher-order coupling is 
essential. One salient ingredient is a proper description of neck formation. Though there 
are many pioneering attempts which relate the large enhancement of the fusion cross sec- 
tion to the neck formation (Jahnke et al., 1982; Krappe et al, 1983; Iwamoto and Harada, 
1987; Iwamoto and Takigawa, 1989; Aguiar et al, 1988), a microscopic description of fu- 
sion reactions in general and neck formation in particular is still at a very primitive stage 
and needs to be further developed. In connection with the former, we wish to note the 
computer simulations for sub-barrier fusion reactions by Bonasera and Kondratyev (1994). 
The effects of neck formation could be formulated in terms of the quantum tunneling in a 
multi-dimensional space (Landowne and Nix, 1981; Kodama et al, 1978; Yu Denisov and 
Royer, 1994; Yu Denisov and Royer, 1995). 

Several existing and presently under construction experimental facilities providing beams 
of short-lived radioactive nuclei present new opportunities to explore the dynamics of the 
fusion reactions below the Coulomb barrier. In such facilities, in addition to testing our 
present understanding of the fusion dynamics in a new setting, one can experimentally 
investigate entirely new facets such as the coupling of resonant states to quantum tunneling 
and the possibility of molecular bond formation. 

One should finally remark that there are many other tunneling phenomena in nuclear 
physics besides heavy-ion fusion reactions. Alpha decay, fission, various rare decays, and 
nuclear structure problems such as the decay of a superdeformed band to a normal band 
will also be effected by coupling to the intrinsic degrees of freedom and insight obtained 
in the study of heavy-ion fusion reactions at subbarrier energies will be a valuable tool to 
understand these phenomena in more detail. 
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FIG. 1. One dimensional potential of Eq. (4) for the 64 Ni+ 64 Ni system for several I values. 
The lowest barrier is for I = (the bare barrier). The middle and top barriers are for I = 100 and 
I = 150 respectively. 

FIG. 2. The effective radius extracted from fusion calculations for the 16 0+ 154 Sm system using 
Eq. (6). The curves correspond to spherical, vibrational and deformed nuclei with quadrupole 
coupling strengths t>2 = 0, 0.13 and 0.26, respectively. From Balantekin et al. (1996). 

FIG. 3. Classical (on the left) and quantum mechanical (on the right) transmission probabilities 
for a one-dimensional potential barrier. 

FIG. 4. Classical (on the left) and quantum mechanical (on the right) transmission probabilities 
for a two-channel coupling. V$ is the height of the one-dimensional potential barrier coupled to 
these channels. 

FIG. 5. Fusion cross section and barrier distribution for the 16 O + 154 Sm system by Leigh et 
al. (1995). 

FIG. 6. Effective one-dimensional potential barriers from (Balantekin et al., 1983). The outer 
turning points are determined from the phenomenological potential of Krappe et al. (1979) to 
fit the peak positions. The distance between the outer and inner turning points is the thickness 
function inverted from the data. The shaded region indicates the error envelope. The short dashed 
line denotes the point Coulomb potential and the long-dashed line denoted the potential of (Krappe 
et al, 1979). 

FIG. 7. Illustration of the geometric interpretation of the sudden approximation. The solid 
curve is the total potential for the 16 0+ 154 Sm system when the projectile is taken to be spherical. 
The dashed (A = —0.327) and dot-dashed (A = +0.613) curves are potentials for the two-level 
approximation for the target with (3 = 0.25. The arrows show the shifts predicted for the barrier 
peaks (Balantekin et al., 1996). 

FIG. 8. A systematic study of subbarrier fusion of 16 projectile with rare earth nuclei, the 
structure of which is described using the Interacting Boson Model (see text). 

FIG. 9. Fusion cross sections from Ackermann et al. (1996) for the systems 58 Ni+ 64 Ni (circles) 
and 64 Ni+ 64 Ni (squares) as a function of the energy normalized to the barrier height. 

FIG. 10. Comparison between measured and calculated barrier distributions for 36 S+ 110 Pd 
system as more phonons are included in the calculation. The data are from Stefanini et al. (1995c) 
and the calculation is from DeWeerd (1996). The dot-dashed, dashed, and solid lines correspond 
to calculations including one-phonon, two-phonon, and three-phonon states, respectively. 

FIG. 11. Fusion cross section (upper panel) and the average angular momenta (lower panel) 
for the 64 Ni+ 96 Zr system. The data is from Stefanini et al. (1992). the theoretical calculation is 
from Hagino et al. (1997a). The two-phonon states of the quadrupole surface vibration of both 
the projectile and the target are taken into account in the coupled channels calculations. The 
dotted line is the result in the absence of channel coupling. The dot-dashed and dashed lines are 
the results when the nuclear potential is expanded up to the first and the second order terms 
in the deformation parameters, respectively. The solid line is the result of the coupled channels 
calculations to all orders, obtained without expanding the nuclear potential. 
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FIG. 12. Predicted behavior of the barrier distributions for fusion reactions on the prolate 
( 194 0s) and the oblate ( 194 Pt) nuclei (Balantekin et ai, 1994a). 



FIG. 13. Experimentally determined fusion cross sections by Bierman et al. (1996a) for prolate 
and oblate nuclei. The solid curve is the simplified coupled channels calculation with the CCDEF 
code. The dashed curve is the result for a one-dimensional barrier ignoring all the couplings. 
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